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to extract the CKM matrix element |y„j| from experimental data on K semileptonic decays. The 
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simulated with the asqtad action (MILC Nf = 2+\ configurations). We discuss the general 
methodology of the calculation, including the use of twisted boundary conditions to get values of 
the momentum transfer close to zero and the different techniques applied for the correlators fits. 
We present initial results for lattice spacings a « 0.12 fm and a ~ 0.09 fm, and several choices of 
the light quark masses. 
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1. Introduction 



The error associated with the lattice determination of the form factor f^^{0) (~ 0.5%) is 
still the dominant uncertainty in the extraction of IV^ I from experimental data on K semileptonic 
decays: |K.v|/+'^(0) = 0.2163(ib0.23%) [|l]]. Improvement in the determination of that form factor 
is thus crucial in order to extract all the information from the available experimental data. 

A precise value of |V„5| is needed to check unitarity in the first row of the CKM matrix. Any 
deviation from unitarity would indicate the existence of beyond the Standard Model physics. But, 
even if unitarity is fulfilled, however, as it is the case with current experimental and theoretical 
inputs, this test can establish very stringent constraints on the scale of the allowed new physics 
(~ 10 TeV) One also could compare the values of \Vus\ as extracted from helicity-allowed 
semileptonic decays and helicity-suppressed leptonic decays in the search for deviations from SM 
predictions. In particular, it is useful to study the ratio 
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where the subscripts pLl and 13 indicate that those quantities are obtained from experimental data 
on leptonic K^2 and semileptonic K^t, decays respectively. The ratio in (|0|) is unity in the SM but 
not in some extensions of the SM, for example, those with a charged Higgs. Again, the error in the 
current value /?^23 = 0.999(7) [|l|] is limited by the precision of lattice-QCD inputs. 

In these proceedings we report on the status of the calculation of the form factor ff'"^ (0) using 
staggered quarks. The goal of this analysis is to show that the staggered formulation can provide 
a determination of this parameter competitive with the state of the art unquenched determinations 
[0] by addressing the main sources of systematic errors and improving in statistics. 



2. Methodology: extracting the form factor directly at 







Semileptonic K decays are parametrized in terms of the form factors /+ and /o in the following 



way 



{7l\V^\K)=fl-{q') 



+ f§^{q')'^^^qr 



(2.1) 



where q = Pk — Pn and V is the appropriate flavour changing vector current. One of the main 
components of our analysis which reduces both systematic and statistical errors is the use of the 
method developed by the HPQCD collaboration to study charm semileptonic decays |^. This 
method is based on the Ward identity relating the matrix element of a vector current to that of 
the corresponding scalar current: q^ {n\'V!l^''\K)Z 



mq){K\S^'"-\K), with S = si, and Z a 
lattice renormalization factor for the vector current. In this work, we use the local scalar density 
of staggered fermions, so the combination (m^ — mq)S requires no renormalization. Using the 



definition of the form factors in Eq. (lA) and this identity, one can extract f^^iq^) at any q^ by 
using 



m.v — mi 
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Figure 1: Structure of the 3 -point functions needed 
to calculate f^^iq^)- Light-quark propagators are 
generated at tsom-ce with random-wall sources. An 
extended sti^ange propagator is generated at T. 
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Figure 2: Deviation of our data from the continuum 
dispersion relation prediction. 



Kinematic constraints demand that /+(0) = /o(0), so this relation can be used to calculate ff^{0). 
One of the main advantages of relation ( |2.2[ ) is that it avoids the use of a renormalization factor to 
obtain the form factor /q. The drawback to this method is that it gives no access to the shape of 
ff^, but in this analysis we are focusing on the extraction of |V„.s|> so we need the normalization of 
the form factor only at a single point. 

Another key ingredient is employing twisted boundary conditions ^ to simulate the rel- 
evant correlations functions directly at = 0. This avoids an interpolation in and thus the 
corresponding systematic uncertainty. In Fig. [T] we plot the general structure of the relevant 3-point 
functions. In order to get = 0, we inject momentum p = dn/L'm either the kaon or the pion. 
For a non-zero pK we chose Bq = 02 = 0, 0i 7^ 0, and for a non-zero pji we chose % = di= 0, 
02/0 (see Fig. [l| for definition of 0o,i,2)- The twisting angles are tuned to produce q^ = using 
two-point correlators fits according to 

We found in a previous test run that the use of random-wall sources greatly reduces the 
statistical errors of the parameters of the two-point and three-point correlators, so we use them 
throughout this analysis. 



3. Simulation details and fitting 

We have completed the generation of correlators with HISQ staggered valence quarks on the 
2 + 1 flavor asqtad staggered MILC ensembles [||] shown in Table |l| at two lattice spacings. We 
average results over four time sources separated by 16 (24) timeslices on the 0.12 fm (0.09 fm) 
ensembles but displaced by a random distance from configuration to configuration to suppress 
autocorrelations. A subset of this data was analyzed in [^. The strange valence mass is tuned to 
its physical value on each ensemble [Q]. The valence light-quark masses are fixed according to the 
relation, = '"iL}"'"^''"''^ ■ The effect of the mixed actions for the sea and valence quark 

nis ' (HISQ) m, ' (asqtad) 

sectors can be analyzed using partially quenched staggered CHPT techniques for the chiral and 
continuum extrapolations. 
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^ a (fm) 


ami /anih 


Nconf 


am™' 


amf 


^ sources 




coarse 0.12 


0.020/0.050 


2052 


0.0491(9) 


0.02806 


4 


5 




0.010/0.050 


2243 


0.0495(9) 


0.01414 


4 


8 




0.005/0.050 


2098 


0.0489(9) 


0.00670 


8 


5 


fine 0.09 


0.0124/0.031 


1996 


0.0337(6) 


0.0080 


4 


5 




0.0062/0.031 


1946 


0.0336(6) 


0.0160 


4 


5 



Table 1: Ensembles and simulation details, am^ is the nominal strange-quark mass in the sea sector, Nsom-ces 
is the number of time sources, and Nj the number of sink-source separations for which we have generated 
data. 



We fit the two-point functions for a pseudoscalar meson P to the expression 

cUPP;t) = I (-ir(zD' {e-^P + e-^^^^^'-'^) , (3.1) 

m=0 

where L, is the temporal size of the lattice. Oscillating terms with (—1)"' do not appear for pions 
with zero momentum. From two-point function fits, we checked whether the continuum dispersion 
relation is satisfied. This is plotted in Fig. ^, which shows very small deviations from the continuum 
prediction (< 0.15%), indicating small discretization effects. 
The functional form for the three-point functions is 

C^p^^'ipMtJ) = f (-l)""(-l)"(^-'U'""Z^Z,f [e-^l(T-t)-EiiT-L,+t) 

m,n=0 

(3.2) 

where the factors Zf are the amplitudes of the two-point functions in (|]l|). The three-point parame- 
ter A°° in (U) is related to the desired form factor f^'^iq^) via f^'^iq^) = ^A^'^ (q^) VlE^ {m^ - 



mi)/ {m\ — m\), where we have used (22) and taken into account some overall factors involved in 
the parametrization of the correlation function. We extract the form factors f^^iq^) using the ex- 
pression above directly from simultaneous fits of the relevant three- and two-point functions. In 
these fits we include several three-point functions with different values of the source-sink separa- 
tion T, with at least one odd T and one even T to be able to get a handle on the contributions from 
the oscillatory states. 

In this analysis it is especially relevant to check for the stability of our fits under the choice of 
fitting parameters and techniques, since we are getting very small statistical errors and we need to be 
sure that these results are not methodology dependent in any way. One of the checks we performed 
is varying the time fitting ranges and number of states included in the fits. Fitting ranges for two- 
point functions are t,m„ — {Lt — tmin) and for the three-point functions — {T — tmin), with Lt the 
temporal size of the lattice and T the source-sink separation — see Fig. The number of states 
included is the same in the regular and oscillating sectors, so Ng^p = A'reguiar states = A'osciiiatory states- 
Fixing Nexp and changing tmin from 3 (5) for 0. 12 fm (0.09 fm) ensembles up to the maximum 
allowed by the source-sink separation, give us a plateau for central values with only small variations 
in errors. Analogously, fixing tmin to our preferred value we do not find any significant variation of 
results for Nexp > 3 — 4. 
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We study as well which combination of T's from the ones we have simulated is optimal. We 
find that the central values are very insensitive to the number of three-point functions included and 
the values of T in the range we are analyzing. Errors and stability are better when 15 < T < 24 and 
including three three-point functions for the 0.12 fm ensembles, and 18 < T < 33 and including 
four three-point functions for the 0.09 fm ensembles. 

Finally, we checked an alternative way of doing the fits, using the iterative superaverage 
method described in [^]. This takes an explicit combination of three-point functions with consec- 
utive values of T and the time slice t which suppresses the contribution from both the first regular 
excited state and the first oscillatory state. Again, results are compatible within one statistical a 
with our preferred fitting method. 




Figure 3: In the left panel, we collect the form factor f^^{Q) obtained from the different ensembles in 
Table [1]. The right panel shows an example of the output from the chiral and continuum extrapolation tests. 
The dotted black line in both plots is the NLO CHPT prediction. The solid black line is the extrapolation 
in the light-quark mass with a physical and in the continuum limit. The green and red lines are also 
extrapolations in the light-quark masses with mf but with a « 0. 12 fm, 0.09 fm, respectively. 



4. Results for /+ ^0) 

In Fig. ^ we collect the results for f^^{0) with statistical errors from our preferred fits as a 
function of {r\mjiY, fixing {2m\ — m\) to the experimental value. We plot results coming from 
three-point functions where the external momentum to obtain = is injected via the K and 
the n (we give an offset in the pion mass to the two points for clarity). The green triangle at the 
far left corresponds to the 0.12 fm ensemble with masses 0.005/0.050, four time sources, and a 
moving pion. When analyzing those data we found it to be more challenging to get stable results, 
so we have decided to double the number of sources and exclude it in our discussions until the full 
new data set is analyzed and other effects, like finite volume corrections, are incorporated to the 
analysis. So, in particular, we do not include it in the fits described below. 

The first remarkable characteristic of our results is that the statistical errors are very small, 
0.1 —0.15%, reaching our goal to be competitive with other determinations. In addition, the results 
coming from three-point functions where the external momentum to get = is injected via the K 
and the 7i agree within the very small statistical errors, as can be seen in the figure. This constitutes 
a very good test of our methodology and quoted errors. 



5 



K — >■ nlv form factors with staggered quarks 



E. Gamiz 



4.1 Chiral and continuum extrapolation 

The form factor /+(0) can be written as a CHPT expansion in the following way: /+(0) = 
1 + /2 + /4 + /e + . . . = 1+ /2 + A/. The AdemoUo-Gatto (AG) theorem, which follows from vector 
current conservation, ensures that /+(0) — )• 1 in the SU (3) limit and, furthermore, that the SU (3) 
breaking effects are second order in {m\ — ni\). This fixes completely in terms of experimental 
quantities. At finite lattice spacing, systematic errors can enter due to, for example, corrections to 



the dispersion relation needed to derive Eq. (2.2). Those and other discretization effects are very 



small though, as can be deduced from our results in Fig. ^ However, since statistical errors are 
at the 0.1 — 0.2% level, we should pin down the other sources of systematic errors as precisely as 
possible. 

Our plan for treating the light-quark mass dependence and the discretization effects in our cal- 



culation is to use two-loop continuum CHPT [10|, supplemented by staggered partially quenched 



CHPT at one-loop. The small variation with a in our data suggests that addressing those effects at 
one loop should be enough for our target precision. 

Since we do not yet have the staggered CHPT expressions, nor have we implemented the two- 
loop continuum CHPT functions, just as an exercise, we can try to fit our data with a much more 



simple fitting ansatz. We take the continuum partially quenched NLO CHPT expression [11] and 
add a general parametrization of NNLO analytic terms and corrections of the form 

(0) = l+/2+C«(^y+ r\{ml-mlf[6^\rym,f + cf\nmKf 

+ di\nm,flog{mll^^)+ct\nm,f + cf 



1. 

'.(4.1) 



We include only correlation functions coming from injecting the momentum in the n for these test 
fits. The result for one of these fits with C„ = = for / = 3, 4 is shown in Fig. ^ We can obtain 
similar good fits with different combinations of terms in (pH]). This must be taken just as a naive 
first try to fit our data and no conclusions should be drawn until we have used staggered CHPT to 
gain more information about the structure of our data. 



5. Conclusions and outlook 

We have completed the generation of the data in Table [I] needed for the calculation of ff^{0). 
Since the time of the conference we have generated data for another coarse ensemble with light- 
quark mass ami = 0.007 to facilitate the chiral extrapolation, which we anticipate is going to be 
our main source of uncertainty. 

The statistical errors in the form factor in all ensembles exceed our expectations, being around 
0.1 —0.15%. We have performed several checks of the robustness of the central fit values and 
errors, by studying the stability with changes in the time range and number of states, the dependence 
on the source-sink separation and number of three-point functions included in the fit, and testing 
alternative methods for fitting the correlation functions. We find it very difficult to make changes 
in the fitting procedure that change the fit results outside the one sigma range. Another very good 
test of our results is the fact that ff^ (0) as extracted from three-point correlation functions with 
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a moving % and a moving K agree with each other. For the final analysis we will redo the fits we 
found to be the optimal, including the correlation functions with both momentum injected in the % 
and the K to further increase the statistics. 

We found very small lattice spacing dependence in our data and the continuum dispersion 
relation is fulfilled at the 0.15% level, but in view of the small statistical error, we plan to study in 
detail the dependence on cP' by using staggered partially quenched CHPT at one-loop. We will also 
investigate the use of two-loop continuum CHPT. 

With all these elements, we expect our calculation to be competitive with the current state of 
the art. 
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